Cluster Persistence: a Discriminating Probe of Soap Froth Dynamics 
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The persistent decay of bubble clusters in coarsening two-dimensional soap froths is measured 
experimentally as a function of cluster volume fraction. Dramatically stronger decay is observed 
in comparison to soap froth models and to measurements and calculations of persistence in other 
systems. The fraction of individual bubbles that contain any persistent area also decays, imply- 
ing significant bubble motion and suggesting that Tl processes play an important role in froth 
persistence. 
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Cellular structures are common in nature PrH, and 
are studied in soap froths, polycrystallites, Potts mod- 
els, biological cell clusters, and magnetic bubble arrays. 
They are composed of cells with varying sizes but identi- 
cal composition, packed tightly against neighbors so that 
cell boundaries meet at vertices. In most cases, this 
structure is dynamic, i.e., part of an ongoing nonequi- 
librium process, and the challenge is to understand the 
role of dynamics in determining the structure. Individual 
cells continually evolve so as to reduce the total interfa- 
cial area in the system, which results in coarsening, or 
growth of the characteristic bubble size. "Topological" 
changes of cell adjacencies occur when neighbors swap 
(Tl process) and when cells are destroyed (T2 process) 
g. Topological transitions are effectively instantaneous 
in soap froths, while bubble areas grow smoothly. This is 
in contrast to polycrystallites or dynamical Potts models 
where topological processes occur at the same slow time 
scales as the evolution of bubble areas. Rapid topological 
transitions in soap froths allow us, in principle, to isolate 
the effects of these processes. 

Soap froths are ideal systems to study cellular pattern 
evolution, since they coarsen at laboratory time scales 
and their time-dependent structure can be directly im- 
aged in two-dimensions (2d) with video microscopy 0. 
In this paper, we study isotropic froths in 2d. We limit 
ourselves to the scaling regime H at late times, when 
the absence of non-universal transients facilitates both 
the analysis and the comparison with model systems. 

A relatively new and general probe of non-equilibrium 
dynamics is the local power-law decay of persistence to- 
wards zero [0, P(t) ~ t~ e , which in froths measures the 
fraction of the system that has always remained within 
the same bubble from time to through to a later time 
t 3> £rj- A simple theoretical argument provides a lower 
bound 6 > 1 : the area of remaining bubbles in a coars- 



ening froth increase linearly in time, implying that the 
number of bubbles N ~ 1/t 011 . Since persistent vol- 
ume cannot increase within a given bubble, then P is 
bounded above by N ~ 1/t asymptotically. Numerical 
studies of 2d cellular systems have found < 1 p[-^0[, 
and are taken in conjunction with this theoretical bound 
to indicate 9 = 1 asymptotically pi. However, exper- 
imental studies of 2d soap froths found 8 > 1 fUJl , in 
clear disagreement, though consistent with the theoreti- 
cal bound of 6 > 1. Evidently, persistence probes inter- 
esting physics in soap froths. 

Most studies of froths use measures specific to the cel- 
lular topology, such as the distribution of the number ver- 
tices per cell. While such measures are important (see, 
e.g., E3]), it is also useful to develop probes of struc- 
ture that are not unique to cellular systems. Results can 
then be compared to corresponding systems that have no 
cellular structure and hence no topological processes, so 
that the distinct effects of topological dynamics can be 
revealed. For example, it is interesting to compare froths 
with non-conserved scalar coarsening since they both ex- 
hibit curvature-driven interface motion while the latter 
has no topological processes. However, the choices for 
such probes are limited, since froth correlations generally 
reflect the cellular structure (see, e.g., |R|) independently 
of the dynamics. 

Two of us have proposed a more general definition of 
persistence that makes use of the soap froth structure 
to create a "non-topological" measure that can be com- 
pared to other systems |jlj,^5| . One first randomly selects 
a fraction </> of all of the bubbles at an initial time to and 
colors them "green." One then asks what fraction of the 
system has always been green between to and t. Essen- 
tially, we first make a virtual phase and then study its 
persistence properties. This allows us to fingerprint the 
dynamics at different length scales, as probed by the size 



'Permanent address: Department of Physics, Dalhousie University, Halifax NS, Canada B3H 3J5 



of clusters of the virtual phase. The limit <\> — + con- 
sists of isolated bubbles whose persistence decays with 
the encroachment of any neighbor, which corresponds to 
the persistent decay previously studied in froths. Phases 
with <j) > have larger clusters, and for sufficiently large 
</> these clusters percolate through the system. Since </> is 
also the volume fraction, this "cluster persistence" can be 
directly compared to persistence in a two-phase system 
at the same volume fraction but lacking an underlying 
cellular structure. 

Our primary results reported in this paper include (i) 
experimental measurements of cluster persistence for all 
volume fractions < <\> < 1 in a coarsening 2d soap froth; 
for all 4> we find a power law decay with exponent 9{<p); 
(ii) these exponents 9{4>) are greater at corresponding 
<p than any other studied system with curvature-driven 
growth, cellular or non-cellular; (iii) the persistence ex- 
ponent is almost linear with </>, for which we provide a 
mean-field argument; and (iv) as (j) — > 1, 9(4>) vanishes 
linearly with a considerably larger slope than predicted 
Jl4| for a mean-field droplet model without cellular struc- 
ture. We interpret these results, particularly point (ii), 
as evidence of the strong role of Tl processes in froths. 

The experimental setup has been reported earlier 
(0,pJJ. We start our experiment by pumping approxi- 
mately 20 000 soap bubbles into a chamber formed by 
two rigid parallel plates (26.7 x 36.8 cm 2 with a 0.16 
cm gap) , making an initially random froth that proceeds 
to coarsen. The froth was well drained before starting 
the experiment. Images of the froth were captured using 
a high resolution CCD camera (1037 x 1344 pixels) ev- 
ery 10 min, sufficient to follow the evolution of individual 
bubbles. [The entire chamber is slightly larger than the 
imaged region. We also excluded strips of bubbles near 
the edge of our image region to test boundary effects: 
none were found.] The scaling regime, as identified by a 
stationary distribution of bubble topologies, set in from 
6 h when there were approximately 6000 bubbles, and 
we continued taking data until 160 h when 200 bubbles 
remained. Our persistence data was taken from to = 12 
h, with 2000 bubbles. Quantitatively consistent results, 
though noisier, are obtained using to — 26 h. 

We analyze the CCD images to identify each pixel with 
a given bubble (see |l7]]). Boundary pixels, on bubble 
boundaries, are randomly assigned to adjacent bubbles. 
For a given subset S of bubbles, we represent the pixels 
in S by a function ps(j,t) = 1, where j G [1, 1393728] 
counts the M pixels, and for pixels outside of that subset 
of bubbles, Ps{j,t) = 0. The cluster persistence of the 
set S is then formally 

M t 

P 5(*) = ^En^(i,0: (1) 

j=l t'=t 

where < P$ < 1- In other words, a pixel counts towards 
the persistence if it has always been in any of the bub- 
bles belonging to S. With this definition, the traditional 



measure of persistence corresponds to the case where S 
consists of a single bubble. More generally we can take 
arbitrary subsets of all of the N(t ) bubbles, where there 
are 2^*°' possible subsets. We sample subsets by ran- 
domly choosing m = <pN(to) bubbles at a time. We fit 
the persistent decay for individual subsets, and then av- 
erage the resulting exponents — binning them with re- 
spect to 4>(<->)- F° r each subset, we use the time-average 
of 4>(S). The resulting persistence decay is characterized 
by a power-law 



P^t) ~ t 
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We expect 6' = 9, since the average bubble area in the 
froth, {A) : is asymptotically proportional to t [||, but 
we measure 9' to eliminate the dependence on the abso- 
lute time origin and to minimize transient effects |10| , |l7f . 
The decay is well fit to a power-law, as illustrated for 
five subsets in the left of fig. [j], and we plot all of the 
resulting exponents in the right of fig. |l| as a function 
of volume fraction, (f>. The errors shown are from the 
statistical variation between subsets at a given volume 
fraction, and can be reduced with more sampling. We 
selected about 1000 multiple-bubble subsets out of the 
22000 possibilities, in addition to the 2000 single-bubble 
subsets. 

In the right of fig. [3 we plot our data for 9{4>) as 
open circles, and include a quadratic fit as a solid curve 
(described below). For comparison, we include data 
from previous studies of corresponding systems with 
curvature-driven growth and volume fraction <f). These 
include the persistence of one phase of the p-state Potts 
model, with cf> = l/p, simulated (|J^] and plotted as tri- 
angles for < <j> < 1/2. Potts systems, particularly when 
p — » oo, have been proposed as soap froth analogues p8|. 
The Potts data fits relatively well on a continuous curve 
(dotted line), interpolating smoothly between p = 2 with 
a two-phase non-cellular structure, 3 < p < oo with a cel- 
lular structure and some bubble coalescence, and p = oo 
with a cellular structure and no bubble coalescence. Also 
shown (open square) is a previous persistence experiment 
on twisted nematic liquid crystals p9[ , an Ising model 
analog with <p — 1/2- Finally, an analytic result for a 
mean- field droplet model, valid for </> « 1, describing 
non-conserved order parameter dynamics with a time- 
dependent applied field is plotted as a dashed line [ JL4J . 
We note that remarkable agreement is seen between the 
Potts data, the experimental work on twisted nematics, 
and the extrapolation of the 0(1 — <fi) exact result for 
coarsening droplets. Neither the nematic nor the droplet 
system have cellular structures, and the relatively good 
agreement between them and the Potts data indicates 
that persistence is not greatly affected by the cellular 
structure per se. Nevertheless, we find froth decay expo- 
nents that are dramatically larger than all these corre- 
sponding systems! 
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FIG. 1. Left: log-log [base 10] plot of cluster persistence decay for various volume fractions cf>. The arrows indicate the fitting 
region used to extract all persistence exponents. Right: persistence decay exponent as a function of volume fraction for our 
soap- froth data (open circles), Potts models with <f> = 1/p (triangles), twisted nematics (open square), and droplet calculation 
(dashed line) . The solid curve is a quadratic fit [see text] , while the dotted curve is a guide to the eye. The arrow indicates the 
4> — value expected from bubble motion, as probed by the decaying number of bubbles with persistent cores [see text and 
fig. 2]. 



Qualitatively, the froth persistence data is strikingly 
linear. The droplet model studied in Ref. JL41 used a 
mean- field treatment, asymptotically correct as (f> — ► 1, 
to obtain the linear variation 6(4>) — 0.48797(1 — 0). This 
motivates a similar mean-field argument for froths. In 
essence, if the local environment around a persistent re- 
gion is entirely random, then the rate of persistent vol- 
ume lost from a subset of bubbles S of volume fraction 4> 
is given simply by the loss of persistent volume from indi- 
vidual bubbles times the chance, l—(f>, that it is being lost 
to a bubble outside of S. That is, P^/P,/, = (1 — <j>)Pq/Pq 
where P$ is the cluster persistence, and Pq is the <f> — > + 
limit of single bubbles. This equation would result in a 
decay exponent 6((f>) — 6(0)(l — 4>). It is a priori un- 
likely that this drastic mean-field approximation, or the 
resultant linearity, holds as an exact relationship. For 
example, in the droplet model |Q appreciable curvature 
corrections, expected to enter at Of(l — (j)) 2 ), must de- 
velop since 6(0) = 1. Indeed, we see below that such 
non-mean-field effects occur in soap froths. However, the 
near linearity of the froth data does imply that correla- 
tion effects are small for cluster persistence in soap froths. 

To fit all of our froth data, we use an expansion in 
powers of 1 — 0, making no assumptions about the be- 
havior near 0=1, and obtain an excellent fit: 6(4>) = 
-O.OOl ±OOO3 + 1.15 ±OO2 (l-0)+O.15 ±OO3 (l-0) 2 . Higher 
powers in (1 — 0) are not indicated, since we find that 
the coefficients are zero within error bars. The quadratic 
term is small but significantly non-zero — this upward 
curving 6(<j)) indicates that the persistence decays, for 
all < 4> < 1, slightly slower than our mean-field ar- 
gument would predict. Evidently cluster persistence re- 



veals positive correlations between persistent regions and 
neighboring cells of the same phase. The variation near 
d> = 1 is linear within errors, as was conjectured in UM 
by analogy with the droplet model. 

We turn now to the question of the strong persistence 
decay in froths relative to other curvature-driven sys- 
tems, for which purpose we consider the standard persis- 
tence exponent, or equivalently the persistence of isolated 
bubbles. The extrapolation to = of our quadratic fit 
gives an improved value of 6(0) — 1.30 ± 0.01 for single 
bubbles p0[ , which significantly disagrees with previous 
predictions of 6(0) = 1. In the scaling state, if a constant 
fraction of remaining bubbles have persistent cores, and 
if these cores have a non-decaying average area, then the 
persistence is proportional to the density, and 6 — 1 . Our 
result of 6 ~ 1.3 implies either a continual erosion of the 
persistent cores in bubbles, a continual decrease in the 
fraction of bubbles with persistent cores, or both. By di- 
rect measurement we find that the average persistent core 
size does not decay at late times, while in contrast, the 
fraction of bubbles with persistent cores steadily decays 
and entirely accounts for the exceptionally strong persis- 
tence decay of individual bubbles in soap froths. This 
strong decay is not shared by solutions of topological or 
droplet mean field models (with or without vertices, re- 
spectively) or by simulations of Potts models. 

We first measure the average persistent area 
(A*(t)/A(t )), where averages are restricted to the N*(t) 
bubbles that contain persistent regions at time t, and we 
have normalized each bubble by its area at tg. As shown 
on the left in fig. 0, it does not change in time after 
50 h, and its distribution approaches a time-independent 



form. This is qualitatively similar to the evolution of di- 
lute droplet models (h|, where (A*(t)/A(t )) w 1 since 
all surviving drops have a persistent core from to. 

We next measure the number of bubbles that contain 
any persistent area, N*(t\, and plot the ratio N* /N vs 
(A(t)) on the right in fig. g. While noisy, there is a clear 
decay, with a best-fit giving N* ~ ^)-i-27±o.02 Tnis 



is a dramatic and unexpected result. It implies a new 
length scale L* ~ i - 64±001 ) characterizing the separation 
of bubbles that have persistent volume, that grows faster 
than the bubble size L ~ t 1 / 2 [|]|3|. It also independently 
implies a persistence decay P * = (A*(£))7V* ~ AT*, indi- 
cated by the arrow in fig. |l| that is consistent with our 
previous result. 
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FIG. 2. Left: plot versus time of {A* (t) / A(to)) , the average of the persistent area within a bubble at time t normalized by 
the entire area of the same bubble at to- The averages are restricted to those bubbles with non-zero persistent regions. The 
horizontal line, at (A*(t)/A(to)) — 0.314 is the average value from 50 — 150 h. The distribution for various times is shown in 
the inset; it approaches a stationary form. Right: fraction of bubbles that contain persistent volume at time t vs. the average 
bubble area. The best fit (solid line) gives a decay exponent —0.27 ± 0.02, where the error bar is approximate. 



The power-law decay of N*(t)/N(t) implies significant 
movement of the center of mass of a bubble, since persis- 
tence within a given surviving bubble cannot be entirely 
lost without the bubble moving off of its original center 
of mass. Bubble motion on the scale of the bubble size 
is therefore directly implicated in the strong persistence 
decay of soap froths. 

It is worthwhile to speculate on the causes of the bub- 
ble motion. We have investigated the decay of persistence 
of individual bubbles, and found numerous discrete de- 
creases in the persistent fraction correlated with signif- 
icant and sudden motion of the bubble center of mass. 
This implicates the correspondingly sudden topological 
processes in bubble motion and hence in strong persis- 
tence decay. T2 processes do not lead to sudden bub- 
ble motion for three-sided bubbles as they vanish, and 
shrinking bubbles with four or five sides loose their ex- 
cess sides through Tl processes before they vanish. As a 
result, Tl processes are the main source of bubble mo- 
tion. Indeed, accurate numerical simulations of planar 
froths have seen that bubble motion extends appreciable 
distances from Tl events [ pi] , in agreement with our ob- 
servations. It had been thought that Tl processes had no 



discernible effect on foam structure due to the unchanged 
growth of bubble size when weak shear was applied p2[ . 
However the growth of bubble size is well predicted by 
mean-field like arguments H that should become more 
accurate with shear due to mixing effects, and so is an 
incomplete probe of structure. In contrast, our results 
on persistence decay show that Tl processes, often not 
included in theoretical work |23j, are important in the 
understanding of soap froth evolution. This is consistent 
with results of numerical modeling that show that the 
width of the distribution of bubble sides, /i2, is depen- 
dent on the specific implementation of Tl processes [ pi| . 

Potts models also have a cellular structure and hence 
topological processes occur throughout their evolution. 
There, within numerical accuracy #'(0) = 1 (§|, and so we 
must have both (A*) and N* /N asymptotically constant 
and hence no significant motion of bubbles. This follows 
from the monotonicity of (A*(t)) and N ~ I /{A) ~ 1/t. 
What leads to this distinction between soap froths and 
p —* oo Potts models? Qualitatively, topological pro- 
cesses occur much more slowly in Potts models as com- 
pared with soap froths 0. However, the relative rate of 
Tl and T2 processes in Potts model evolution has neither 



been investigated nor compared with soap froths, where 
the ratio of T1::T2 is approximately 3 :: 2 pi[ . As yet, 
it is impossible to distinguish the role of the frequency of 
topological T\ processes from their speed in the strong 
decay of persistence. 

In summary, from 2d experiments of coarsening soap 
froths, we have measured persistence decay exponents as 
a function of volume fraction. We find that soap froths 
exhibit much faster persistence decay at all volume frac- 
tions compared to other systems, including existing mod- 
els of soap froth evolution in which persistence has been 
measured. Over time, a decreasing fraction of surviving 
bubbles contain any persistent area at all. This implies 
significant motion of bubbles that is strongly correlated 
with topological processes in the froth. 

While direct study of T\ processes is clearly desirable, 
we have demonstrated that the decay of persistence is a 
sensitive probe of the effect of T\ processes. More fun- 
damentally, persistence decay also allows the comparison 
of different cellular systems, for instance between planar 
and bulk soap froths. In light of the possible prolifera- 
tion of T\ processes in the late-time regime of bulk 3d 
froths |25| , it would be interesting to analyze persistence 
decay in those systems. Such an exercise is now possible, 
in principle, with the development of spatially resolved 
MM techniques |||. 
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